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ABSTRACT 

A radiative transfer method appropriate for use in simple cli- 
mate models and 3-D global climate models has been devel- 
oped. It is fully interactive with climate changes, such as 
in the temperature-pressure profile, cloud distribution and 
atmospheric composition, and it is accurate throughout the 
troposphere and stratosphere. The vertical inhomogeneity 
of the atmosphere is accounted for by assuming a correla- 
tion of gaseous k-dlstributions of different pressures and 
temperatures. Line-by-line calculations are made to de- 
monstrate that the method is remarkably accurate. The 
method is then used in a 1-D radiative-convective climate 
model to stu<fy the effect of cirrus clouds on surface tem- 
perature. It is shown that an increase in cirrus cloud cover 
can cause a significant warming of the troposphere and the 
earth ^s surl^ce, by the mechanism of an enhanced green- 
house effect. The dependence of this phenomenon on cloud 
optical thickness, altitude and latitude is investigated. 

Introduction. Radiation calculations in climate models should be 
designed to be properly responsive to climate changes, such as in the at- 
mospheric temperature and in the cloud, aerosol and gaseous atmospher- 
ic composition. The major difficulties are (1) accurately integrating over 
complex and sometimes overlapping absorption bands, which greatly 
change their characteristics with height, and (2) accurately including 
the effects of multiple scattering. 

The basis of the approach we take for both solar and terrestrial 
radiation, which we call the correlated k-distrlbution method, is a gen- 
eralization of the k-distrlbution method used by Lads and Hansen (1974). 
Jn this generalization the vertical inhomogeneity of the atmosphere is 
approximately accoimted for by assuming a simple correlation of k-dls- 
tributions at different temperatures and pressures. By means of Hine- 
by-llne’ calculations we demonstrate that the method is remarkably ac- 
curate even for the notoriously difficult 9. Gym ozone band. 

There are several advantages to this radiative treatment, in addi- 
tion to the fact that it includes all significant atmospheric constituents. 
Realistic spectral properties of clouds and aerosols are employed, based 
on Mle scattering computations with the best available composition and 
refractive index information; thus the treatments in the solar and thermal 
regions are self-consistent and we avoid crude assumptions such as black 


309 


or *h^-black* clouds. The accurate treatment employed for solar zen- 
ith angle effects is necessary for assuring proper latitude variations of 
radiative heating. 

Correlated k-distrihution method. The k-distribution f(k) for a 
given gas and frequency interval is the probability density function such 
that f(k) dk is the fraction of the frequency interval for which the ab- 
sorption coefficient is between k and k + dk. The basic idea of grouping 
frequency intervals of gaseous spectra according to absorption coeffici- 
ent strengths goes back at least to Ambartsumian (1936) who used it in 
estimating the influence of absorption lines in stellar atmospheres. 

Lacis and Hansen (1974) used k-distributions to include the effects of 
multiple scattering; they numerically derived a mean atmospheric k- 
distribution for water vapor absorption and scaled the gas amount as a 
function of pressure and temperature to crudely approximate the effect 
of atmospheric inhomogeneity. 

Our generalization of the k-distribution method provides a more 
accurate treatment of the vertical inhomogeneity of the atmosphere; it 
is based on the assumption that the k-dlstributions at all altitudes are 
simply correlated in frequency space, i. e. , the strongest absorption 
occurs at the same frequencies at all altitudes and similarly for the 
weakest absorption. This maintenance of the same relative rank of ab- 
sorption coefficients over the entire pressure range of the atmosphere 
is rigorously correct for a single spectral line with a fixed center and 
for a uniform Elsasser (1942) band model. Fig. 1 is a schematic illus- 
tration of the k-distribution method. 

For a real gaseous spectrum, we expect some blurring of the as- 
sumed correlation as a result of partial overlapping of lines with differ- 
ent line-widths and strengths and the fact that ^e temperature depend- 
ence of the line-strengths is not the same for all lines. To demonstrate 
that this has little impact on the overall accuracy of our method we show 
in Fig. 2 the thermal cooling rates computed with the correlated k-dls- 
tribution method and with line-by-line calculations for the 9. G^m ozone 
band. 

The computations for Fig. 2 were made for standard atmospheric 
temperature and ozone distributions with a 1 km vertical resolution 
(Lacis and Wang, 1979). The solid line was obtained by means of line- 
by-line computations using the line coefficients tabulated by McClatchey 
et al. (1973) with approximately 2 x 10® frequency intervals. The dotted 
curve was obtained by means of the correlated k-distribution approach 
with 5 probability intervals, i. e. , 5 discrete values of k, each with an 
appropriate pressure-temperature dependence. Since the computing 
time is proportional to the number of frequency intervals or k values, 
the relative speed of the correlated k-distribution method is apparent. 

The k-distribution is formally related to the transmission func- 
tion by 

oo 

T(u) = ^ J f(k)e“*“dk, (1) 
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where u is the gas amount. Because of Its convenient mathematical pro- 
perties, we choose the Malkmus (1967) model to represent T(u), 

T(u) = exp { - B [(1 + 4 Su/tiB)^ - 1]} (2) 

where the two parameters, S and B, are the effective line strength and 
line width. Since f(k) is the Laplace transform of T(u) we obtain the k- 
distribution directly in terms of the Malkmus model parameters, 

f(k) = ^ k"^^^(SB)^ exp[^ B(2 - k/S - S/k) ] . (3) 

The Malkmus model k-distribution is analytically integratable to yield 
the cumulative probability 

g(k) = f f(k’)dk». (4) 

This makes it convenient to rank intervals of Ag according to their rela- 
tive absorption coefficient value. This procedure for ranking Ag retains 
an implicit frequency correlation with altitude and provides the basis for 
accurate computation of multiple scattering effects in an inhomogeneous 
atmosphere. 

For a given gas and frequency interval, the k-distribution can be 
obtained in several different ways (Lacis and Wang, 1979). We use the 
tabulated line coefficients for H 2 O, C02t O 3 , CH 4 and N 2 O compiled by 
McClatchey et al. (1973) to compute line-by-line absorption over homo- 
geneous paths for pressures and temperatures over the range of values 
encountered in the stratosphere and troposphere. For each gas we ob- 
tain a table of S and B parameters by least square fitting of the Malkmus 
model to the line-by-line computations. This yields band model para- 
meters which fit the line-by-line results within 1 or 2 per cent over a 
pressure-temperature range from the ground to a height of "^70 km, 
where the assumption of local thermodynamic equilibrium begins to de- 
teriorate. Note that the numerical approach for extracting band model 
parameters is not limited by the changing line shape with altitude; we 
can use the appropriate Voigt profile for each altitude and obtain a cor- 
responding k-distribution in terms of its effective line strength and line 
width parameters. 

Climatic effect of cirrus clouds^ It is well known that clouds 
have two major competing effects on global climate: ( 1 ) they reflect solar 
radiation, which tends to decrease the surface temperature, and ( 2 ) they 
blanket thermal radiation from the earth ^s surfeice, r eradiating at a low- 
er temperature and thus tending to warm the surface temperature. The 
opinion of clionatologists seems to be that a general Increase in cloudi- 
ness would decrease the surface temperature, for example, according to 
Schneider (1972): . .the effect of a sustained increase in the average 

amount of cloud cover of the earth would be a decrease in the global -av- 
erage surface temperature, provided that the cloud top height and cloud 
albedo remain unchanged. 

It is recognized that cirrus clouds can present an exception to the 
expected cooling effect of clouds. For example, Fig. 3 (from Manabe 
and Strlckler, 1964) shows that cirrus clouds at high altitudes can have a 
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Fig. 3. Thermal equilibrium of var- 
ious atmospheres with clouds. The 
height of overcast clouds used for 
each computation is shown by , H 2 
and H 3 for hi^ clouds, M and L for 
middle and low clouds. The equUi- 
brium curve of the clear atmosphere 
is shown by a thick dashed line. (After 
Manabe and Strlckler, 1964), 



OPTICAL THICKNESS (X=0.55^m) 

Fig. 5. Albedo and blackness of cir- 
rus clouds as a function of optical 
thickness, computed for ice parti- 
cles of 25 ym effective radius. The 
^blackness' is the sum of the cloud 
absorbtivity and reflectivity for 
blackbody radiation of 250 K; the 
reflectivity contributes at most 
0.05 to the blackness. 


Fig. 4. Measured cirrus cloud 
infrared emissivities and the 
"critical blackness" curve from 
Manabe and Strlckler (1964). 
(After Cox, 1971). 



Fig. 6 . Effect of cirrus clouds on 
surface temperature (in degrees 
K) as a function of the cirrus op- 
tical thickness and cloudtop 
height. The cirrus are assumed 
to cover 10% of the sky. Calcula- 
tions are for a global average 
temperature profile, which is 
also similar to the mean profile 
at 30-4(f latitude. 
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heating effect, though the magnitude is small compared to the cooling by 
other clouds. According to Manabe and Striclder (1964): ^’If the height of 
cirrus clouds Is greater than 9 km and its blackness for infrared radia- 
tion is larger than 50 percent, cirrus has a heating effect on the temper- 
ature of the earth’s surfece." The cirrus Infrared emissivitles measur- 
ed by Cox (1971) also suggest that a warming effect by cirrus is small 
and rather unusual, cf. , Fig. 4 his measurements are Interpreted as 
indicating that in the tropics only cirrus at altitudes 250-400 mb cause 
heating, and at midlatitudes all cirrus cause cooling. 

However, the definition for blackness used Cox (1971) differs 
from that of Manabe and Strickler (1964) which thus leaves the conclu- 
sions of Cox open to question. The nature of the calculations of Manabe 
and Strickler was to use a fixed cirrus albedo of 21% and to look at how 
the surface heating or cooling would change for different assumed infra- 
red emissivitles. In general the calculations of Manabe and Wetherald 
are extremely valuable,but in order to Investigate the effect of cirrus on 
the temperature profile it is essential that the solar and thermal spec- 
tral regions be treated in a self-consistent fashion using spectrally de- 
pendent optical properties of water ice. The method described above 
permits tills computation to be done with reliable accuracy. 

Results for the solar albedo and thermal ’blackness' are shown in 
Fig. 5 as a function of optical depth. It is apparent that the thermal 
blackness (defined to be 1 minus the total transmission) is much more 
rapidly saturated than is the solar albedo. Typical cirrus cloud emls- 
sivities of '^0.4 to 0.6 measured by Cox (1971) correspond to optical 
depths of ^0. 5 to 1. 0 and should have corresponding albedos of^^5% to 
^10% rather than the widely used value of 21% estimated by Haurwitz 
(1948), As shown in Fig. 6, cirrus clouds have a substantial warming 
effect on the surface temperature. It is only when the optical thickness 
becomes very large and the clouds appear at low altitudes that a cooling 
effect is experienced. 

The most effective optical thickness for the heating (greenhouse) 
effect is approximately x = 2, for which the thermal blackness is nearly 
saturated while the solar albedo is still rather small. The warming ef- 
fect is also sensitive to cloudtop height; thus high clouds, which reradi- 
ate the absorbed surface radiation at a lower effective temperature, pro- 
vide greater warming than low clouds. 

Discussion. These results, which were computed for mean glo- 
bal atmospheric conditions, show a substantial warming effect due to 
cirrus clouds. The same qualitative conclusion holds for tropical and 
polar latitudes. The magnitude of this greenhouse effect is stiff iciently 
large to indicate the need for a careful reappraisal of the role of cirrus 


clouds in climate. 
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